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POTENCY  PROBIT  ANALYSIS 

62507  ! 

This  program  obtains  the  potency  (ratio  of  equally  effective  doses)  and  its 
95%  confidence  limits  for  upwards  of  five  parallel  dose-response  lines  through 
a regression  of  log-dose  on  probit  response  employing  the  maximum  likelihood 
procedure  described  by  D.  J.  Finney  in  Probit  Analysis  (Cambridge  University  Press, 
1952). 

Input  consists  of  a parameter  card  for  the  experiment  followed  by  up  to 
five  subsets  of  data.  Each  subset  of  data  is  preceded  by  a subparameter  card 
and  each  subset  of  data  may  contain  from  two  to  twenty  data  cards  (doses). 

The  parameter  card  for  the  experiment  consists  of  three  4-digit  numbers  which 
identify  the  "experiment",  a dummy  number  identifying  the  first  subset  of  data, 
and  the  number  of  subsets  of  data  in  the  experiment.  In  addition  to  this,  the 
first  parameter  card  also  contains  two  Chi  square  values  (for  testing  for 
homogeneity  of  the  data  and  parallelism  of  the  regression  lines  respectively), 
the  value  of  Student’s  "t"  for  setting  confidence  limits  for  heterogeneous  data, 
the  critical  value  of  for  testing  significance  of  the  regression,  and  a 
"Log  Factor"  for  making  all  doses  greater  than  unity  before  transforming  dose  to 
log^Cdose).  Each  subset  of  data  that  follows  (up  to  five)  must;  be  preceded  by 
a sub-parameter  card. 

Each  sub-parameter  card  contains  three  4-digit  numbers  which  identify  the 
experiment,  the  subset  of  data  to  which  the  sub-parameter  card  belongs,  the 
number  of  data  cards  (doses)  composing  the  subset,  the  number  of  animals  that  were 
observed,  and  the  number  that  responded  in  each  untreated  control  (check) . The 
use  of  sub-parameter  cards  permits  a different  check  for  each  subset  of  data. 

Each  data  card  contains  two  4-digit  numbers  identifying  the  experiment  and 
the  subset  to  which  the  data  belong  (for  ease  of  sorting  and  ordering  of  the 
data  cards),  followed  by  three  6-digit  fields  which  contain  the  dose  used,  the 
number  of  animals  treated,  and  the  number  that  responded. 
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The  program  computes  the  percent  response  in  the  untreated  control,  with 
provision  for  zero  (0)  observations  and  /or  zero  response  in  the  untreated  control. 
The  program,  as  it  reads  each  data  card,  converts  the  dose  to  log^g(dose)  and  the 
numbers  observed  and  responded  to  a proportion,  adjusting  each  proportion  for 
the  untreated  control  by  use  of  Abbott's  Formula.  The  adjusted  (net)  responses 
are  then  transformed  to  probits  by  use  of  a polynominal  approximation.  (Hastings 
1955,  Approximations  for  Digital  Computers).  The  program  then  types  out  the 
numbers  identifying  the  experiment,  the  subset  to  which  the  observation  belongs, 
the  dose  administered,  the  net  response  (as  a percent)  and  the  proportion  that 
responded  in  the  untreated  control.  After  typing  out  this  information  the 
weighted  regression  of  log-dose  on  probit  response  is  computed  using  an  iterative 
procedure  (maximum  likelihood)  with  the  restriction  that  the  dose-response  lines 
for  all  subsets  are  parallel. 

Iteration  continues  until  the  last  and  next  to  last  slope  values  (B's) 
agree  within  ±.00005  or  until  20  cycles  have  been  completed.  After  cycling  until 
convergence  the  program  types  out  the  numbers  identifying  the  experiment  and 
subset,  the  number  of  doses  in  each  subset,  the  weighted  and  corrected  (for 
weighted  mean)  sums  of  squares  of  log-dose  (X),  probit-response  (Y),  cross  products 
(XY),  and  the  sum  of  squares  attributable  to  regression  for  each  subset  (dose- 
response  line).  The  program  continues,  again  typing  out  the  number  identifying  the 
experiment,  the  number  of  doses  in  each  subset,  the  weighted  means  of  log-dose  and 
probit-response,  the  sum  of  the  weighting  coefficients  for  each  subset  and  the 
intercept  of  each  weighted  probit  regression  line. 

Following  this  the  program  next  presents  a summary  of  the  experiment  in  an 
analysis  of  linear  regression  format,  which  includes  the  degrees  of  freedom 
(computed),  sum  of  squares  (total,  regression,  deviations  from  regression,  and 
parallelism),  and  pertinent  parameters  (F  ratio  needed,  F ratio  obtained,  Chi 
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square  values  used  in  the  tests  for  homogeneity  and  parallelism).  Parallelism  is 
the  difference  between  the  sum  of  squares  attributable  to  regression  resulting 
from  fitting  ixxdividual  regression  lines  to  each  subset  of  data  and  from  fitting 
a single  regression  line  (the  same  slope)  to  all  subsets  of  data.  This  definition 
of  parallelism,  because  it  is  an  over  simplification,  is  not  entirely  correct 
as  stated,  and  the  reader  interested  in  the  exact  procedure  is  referred  to  either 
the  accompanying  flow  sheet  or  to  D.  J.  Finney's  Probit  Analysis, 

Following  the  presentation  of  the  analysis  of  linear  regression,  the  program 
then  test  for  homogeneity  of  the  data  before  typing  out  the  number  of  iterations 
completed  (first  line  obtained  from  a weighted  regression  of  log-dose  on  probit- 
response  using  the  observed  proportions  that  responded  in  the  weighting  coef- 
ficients and  is  not  counted  in  determining  the  number  of  iterations),  g (the 
precision  of  estimating  the  slope  which  enters  in  setting  confidence  limits  about 
the  potencies),  the  critical  value  of  used  in  computing  the  confidence  limits 
about  the  potencies,  the  standard  error,  the  standard  error  of  the  slope,  and 
the  slope  of  the  regression  lines  (single  slope  for  all  regression  lines).  If  the 
data  are  heterogeneous,  is  set  equal  to  Student's  with  degrees  of  freedom 

associated  with  the  deviations  mean  square  (total  degrees  of  freedom  minus  the 
number  of  lines,  where  the  total  degrees  of  freedom  is  the  sum  of  the  number  of 
observations  in  the  experiment  minus  one  for  each  dose-response  line  within  the 
experiment),  and  the  square  root  of  the  deviations  mean  square  is  set  equal  to  the 
standard  error.  If  the  data  are  homogeneous  t^  “ 1.96  and  standard  error  =1,00 
are  used  in  computing  the  95%  confidence  limits  of  the  potencies.  The  standard 
error  of  the  slope  (STDEB)  is  invariable  °J ( standard  error) ^/SS  . 

The  program  then  test  for  parallelism.  If  the  lines  are  not  parallel 
LIKES  NOT  PARALLEL  is  typed  out,  and  the  program  proceeds  without  pause  to  read 
the  parameter  card  for  the  next  experiment,  bypassing  the  remaining  calculations 
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in  the  program,  if  the  lines  are  parallel  the  program  continues  with  the  remaining 
calculations,  testing  next  for  significance  of  the  regression. 

If  the  regression  is  not  significant,  as  tested  by  the  critical  value  of  F^,., 
NON  SIG  REGRESSION  is  typed  out  and  the  program  proceeds  without  pause  to  read 
the  parameter  card  for  the  next  experiment.  If  the  regression  is  significant  the 
program  continues,  computing  and  typing  out  the  LD-3Q,  LB-50,  LD-70  and  LD-90  for 
each  regression  line  (subset  of  data).  The  program  then  proceeds  to  compute  and 
type  out  potencies  (ratios  of  LB-50  of  first  line  divided  by  LB- 50  of  each  succeed- 
ing dose-response  line)  and  their  95%  confidence  limits.  The  program  then  proceeds 
without  pause  to  read  the  parameter  card  for  the  next  experiment,  which  initiates 
another  series  of  calculations  for  fitting  parallel  dose-response  lines  using  the 
probit  analysis  procedure  outlined  by  D.  J.  Finney  in  Probit  Analysis. 

The  weighting  coefficients  and  the  procedures  used  in  the  program  are  similar 
to  those  used  in  a Probit  Analysis  program  written  by  R.  J.  Baum,  Clyde  Givens  and 
Gary  Bearden  (Biometrical  Services,  USDA,  Beltsville,  Md„)  for  fitting  simple, 
individual,  dose-response  lines  using  the  Maximum  Likelihood  procedure  described  by 
D.  J.  Finney  (ibid).  The  reader  seeking  greater  detail  on  these  points  is  referred 
to  this  program.  Running  time  with  this  program  is  approximately  four  minutes  with 
typewriter  output.  If  it  is  anticipated  that  more  than  an  occasional  set  of  data 
will  be  processed  by  this  program  it  is  recommended  that  a card-punch  output  be 
used.  This  change  may  be  accomplished  by  simply  changing  the  PRINT  statements  to 
PUNCH  and  recompiling  the  program. 

The  advantages  of  this  program  are  (1)  a neatly  labeled  and  well  organized 
print-out  of  sufficient  information  to  permit  plotting  the  original  and  computed 
dose-response  lines,  and  to  spot  the  subset  (s)  of  data  which  resulted  in  either 
biological  or  statistical  invalidity  if  such  occurs,  (2)  the  inclusion  of  numerous 
checks  which  permits  the  analysis  to  be  completed  only  on  data  which  are  both 
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statistically  and  biologically  valid,  and  (3)  the  use  of  procedures  which  permit 
an  indefinite  number  of  analyses  to  be  performed  without  additional  instructions 

or  interference  from  the  machine  operator® 
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POTENCY  PROBIT  ANALYSIS  MAXIMUM  LIKELIHOOD,  1620  40K,  MULTIPLE  6 
FORMAT  FORTRAN,  AUTO  DIVISION.  THIS  PROGRAM  OBTAINS  THE  WEIGHTED 
LINEAR  REGRESSION  OF  PROBIT  RESPONSE  ON  LOG  DOSE  BY  THE  MAXIMUM 
LIKELIHOOD  PROCEDURE  DESCRIBED  BY  D.J.  FINNEY  IN  PROBIT  ANALYSIS 
(CAMBRIDGE  UNIVERSITY  PRESS).  POTENCY  IS  THE  RATIO  OF  THE  FIRST 
LD-50  DIVIDED  BY  THE  LD— 50  OF  EACH  SUCCEEDING  SET  OF  DATA  IN 
SAME  ORDER  AS  DATA  ARE  READ  INTO  COMPUTER.  MAXIMUM  OF  5 SETS 
OF  20  DOSES  EACH. 

AUTHORS  R.J,  DAUM  AND  CLYDE  GIVENS,  BIOMETRICAL  SERVICES,  USDA, 
BELTSVILLE  MARYLAND,  MAY  1963. 

DIMENSION  C(5),K(5),VN(20,5),X(20,5),P(20,5),WN(20,5),A(5),XB(5) 
DIMENSION  YB(5) ,SX(5),SY(5) ,SXY(5) ,SWN(5),Y(20,5) 

49  FORMAT  ( 2 14, 3F6- 0) 

50  FORMAT  (43H  POTENCY  PROBIT  ANALYSIS  MAXIMUM  LIKELIHOOD,///) 

52  FORMAT (22H  ! DENT  SET  DOSE , 1 OX , 1 2HNET  RESPONSE , 7X , 5HCHECK) 

53  FORMAT  (12HSET  NO. DOSES , 6X, 3HSSX , 1 3X , 3HSSY , 1 3X , 4HSSXY , 1 2X , 5HSSREG ) 

51  FORMAT  (/) 

54  FORMAT ( 1 2HSET  NO. DOSES , 6X , 4HXBAR , 1 2X4HYBAR , 1 2X , 3HSNW, 9X , 9H I NTERCEPT ) 

55  FORMAT  ( / 3 OH  ANALYSIS  OF  LINEAR  REGRESSION/) 

56  FORMAT  (20H  SOURCE  VARIATION  DF, 1 1 X , 2HSS , 1 4X , 2HMS , ) 


57  FORMAT  ( 1 8 H TOTAL 
58  FORMAT  (18H  REGRESSION 

59  FORMAT  (18H  DEV  REGRESSION 

60  FORMAT  (3I4,5F6.0) 

62  FORMAT  (2I5.3F16.6) 

63  FORMAT  ( 2 I 5,4F1 6. 6) 

64  FORMAT  ( i 3 , 5F1 2.6) 

70  FORMAT  (18H  PARALLELISM 


14, F16.6, 16X,F13.3,6H  F05) 

I4,2F16.6,F13.3,6H  FCAL) 
I4,2F16.6,F13.3,8H  CHI  SQ) 


,14, 2FT 6 . 6 , FI 3 . 3,8H  CHI  SQ) 


71  FORMAT  (3H  IT,8X, 1HG,9X,3HTO5,8X,4HSTDE,8X,5HSTDEB,7X,5HSL0PE) 

72  FORMAT  (18HN0N  SIG  REGRESSION) 

73  FORMAT  (18HLINES  NOT  PARALLEL) 

/4F0RMAT ( 1 OH  I DENT  SET , 9X , 4HLD30, 1 2X , 4HLD50, 1 2X , 4HLD7 0, 1 2X , 4HLD9 0) 
75F0RMAT(10HIDENT  SET , 4X , 1 1 HUPPER  L IM I T , 7X , 7HP0TENCY , 7X , 1 1 HLOWER  LIMIT) 
D 1 =7 • 0523 0/8 E- 02 
D 2=4. 2282  01  2E-02 
03=9.2/ 05272E-03 
D4=1 .5201430E-04 
D5=2. 7656720E-04 
D6=4. 30638  00E-05 

1 READ  60, IDE, IDS,N,CH,CH2,T5,F5,FAL 
PRINT  50 
PRINT  52 
I T =—  1 

DO  11  J=1 , N 

READ  60, IDE, IDS,M,B2,B 
PRINT  51 

I F(B2)99, 1 01 ,102 

101  C(J)  = 0. 

GO  TO  103 

102  C ( J )=B/B2 

1 03  G=1 ./(I .-C( J) ) 

K( J )=M 

DO  11  1=1, M 

READ  49, IDE, I DS , Z , VN ( I , J ) , Q 
Q=( (Q/VN( I ,J))-C(J))*G 
B=Q*  1 00. 


10778 
T i08  98 


Til  34 
11158 
Til  66 
YU  234 
11253 
T 1 326 
11362 
11434 
T 1 506 
T 1!  638 
T1722 
T2  022 
T2150 
12246 
12254 
12278 
12374 
12446 
1247  0 
1 24 94 
T2550 
12574 
12598 
12622 
1 2646 
12670 
12694 
12706 
T 2754 
128  02 
12850 
12898 
12946 
12994 
I.3018 
13042 
13066 
13114 
13126 
1321  0 
13246 
13354 
T3498 
1.3658 
13618 
1 3666 
13690 
13714 
13738 
13750 
13834 
13954 


X( i , J)=.4342945*LOG(Z*FAL) 

P( ! s J )™Q 

PRINT  62, IDE, !DS,Z,B,C(J) 

IF  (0)99,3,4 

3 Q”, 0001 

GO  TO  8 

4 SF  (Q-1.)6,5,99 

5 Q~oS999 

6 !F  (Q—  .'5)8 ,8,7 

7 Q-1.-Q 

8 E=SQRT ( LOG ( 1 . / ( Q*Q ) ) ) 

B=2.515517+E*(E*.O1O328+.802853) 

E“E~(B/( 1 .+(E*(E*(E*. 001308+. 189269)+1 .432788)))) 
Z~e  3989421 5*EXP(-E*E*. 5) 

WN( 1 ,J)=(VN( I , J)*Z*Z)/(Q+((1.-Q)+C(J)*G)) 
IF(P(I,J)-. 5)9,9, 10 

9 Y( I , J)=5.~E 
GO  TO  11 

99  IT=1 

10  Y(!sJ)-5.+E 

11  CONTINUE 
PRINT  51 
B=  0. 

SF  (11)12,1,1 

12  SSY-0, 

SSX-O. 

SSXY-O. 

SSREG=0. 

SSNW=0„ 

B2"B 

DO  15  J=1 , N 
M=K(J) 

XB(  J)  = 0. 

YB ( J )=0. 

SX(J)  = 0. 

SY(  J)  = 0. 

SXY  ( J )--0„ 

A1«0. 

A2=0O 

A3=0. 

SWN(  J ) = 0. 

DO  13  1=1, M 
B~WN( I , J ) 

A1=A1 +B 

A2=A2+B*X( I , J) 

13  A|=A3+B*Y( I , J) 

XB( J)=A2/A1 
YB(J)=A3/A1 
SWN( J)=A1 

A1  = 0. 

A2=0« 

A3"0« 

DO  14  1=1, M 
B=WN( S , J) 

Q=X( I , J)-XB( J) 

E=Y( I , J)-YB( J) 
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T4  0/4 

14134 

14194 

14290 

14338 

14386 

14434 

14494 

J4530 

1 4566 

14602 

14698 

14/34 

T4/70 

T48  06 

T48/4 

14942 

15022 

15102 

15 1 14 

15186 

T5234 

J5366 

15378 

1551  0 

15558 

15626 

15662 

15670 

T5694 

T5702 

J5858 

15954 

16022 

16058 

1 6 066 

T61  02 

16234 

1641 4 

16786 

T6794 

T 68 1 8 

1684-2 

T6878 

T6890 

17010 

17  082 

17298 

17322 

17358 

17394 

17430 

T7466 

17490 

17514 

T7538 


A 1 — A 1 +B*Q*Q 
A2«A2+B*E*E 

14  A3»A3+B*E*Q 
SX(J)=A1 
SY(J)=A2 
SXY(J)=A3 
SSNW=SSNW+SWN(J) 

S$X=SSX+A1 

SSY=SSY+A2 

SSXY-SSXY+A3 

15  SSREG=SSREG+A3*A3/A1 
B=SSXY/SSX 
SREG=B*SSXY 

i T= I T+1 

I F{ I T-20) 222, 122,122 
222  ! F (B2-B)23,25,24 

23  IF  (B2-B+. 00005)122, 25, 25 

24  IF  ( B2-B-, 00005)25 ,122,122 
122  DO  22  J=1 , N 

G=1 ./( 1 .-C(J) ) 

H=K( J ) 

A(J)=YB(J)-B*XB(J) 

DO  22  1=1, M 
Q=A(J)+B*X( I ,J) 

E=(Q-5«)/1 .4142136 
IF  (Q— 5- ) 1 6, 1 7> 18 

16  E=-E 

GO  TO  18 

17  E=. 5 

GO  TO  21 

1 8 E= E* ( E*  ( E*  ( E*  ( E*  ( E*D6+D5 ) +D4 ) +D 3 ) +D 2 ) +D 1 ) 

E=(1.-1./(1.+E)**16.)*.5 
IF  (Q-5J19, 21,20 

19  E= , 5~E 
GO  TO  21 

20  E= . 5+E 

21  Z=.39894215*EXP(-(Q-5.)*(Q-5.)*.5) 

Y( I , J )=Q+( (P( I ,J)-E)/Z) 

22  WN( I ,J)=(VN( I ,J)*Z*Z)/((1 .-E)*(E+C(J)*G)) 
GO  TO  12 

25  PRINT  53 
LTDF=0 

SSPAR=SSREG-SREG 
DO  26  J=1,N 

A3=SXY(J)*SXY(J)/SX(J) 

LTDF=LTDF-i-K(  J )-1 

26  PRINT  63,J,K(J),SX(J),SY(J),SXY(J),A3 
LRDF=1 

LPDF=N-1 

LDDF=LTDF~N 

A2=LDDF 

A3=LPDF 

C(4)=6. 28 1 6 

C(3)=5.5244 

C( 2)=5 • 

C(1)=4,4756 


17562 

SPAR=SSPAR/A3 

17598 

SSDEV=SSY-SSREG 

17634 

SDEV=SSDEV/A2 

T'767  9 

FCAL=SREG/SDEV 

T7706 

PRINT  51 

T7730 

PRINT  54 

1775% 

DO  27  J=1 , N 

T7766 

27  PRINT  63, J,K(J) ,XB(J) ,YB( J) 

1 b 00b 

PRINT  55 

18  030 

PRINT  56 

18  054 

PRINT  57, LTDF, SSY , F5 

T8 1 02 

PRINT  58 , LRDF, SREG , SREG , FCAL 

IS  162 

PRINT  59, EDDF, SSDEV, SDEV, CH 

18222 

PRINT  70, LPDF, SSPAR, SPAR, CH2 

1828  2 

1 F( SSDEV-CH) 28 ,29,29 

18350 

28 

Z=1. 

T8374 

T 5= 1 *96 

18398 

GO  TO  30 

T8406 

29 

Z=SDEV 

18430 

30 

SSXY=SQRT ( Z/SSX ) 

T84?8 

G=T  5*T  5*Z / S REG 

18538 

Z-SQRT(Z) 

18562 

PRINT  51 

? 8 58  6 

PRINT  71 

1861  0 

PRINT  64, 1 T , G , T5 , Z , SSXY , B 

18  694 

PRINT  51 

18718 

IF  (G-. 05)31, 32, 32 

T8786 

31 

G=Q* 

.81 

°8 

"o 


8302 
831  0 


13034 
T 9 046 
T3058 
J9262 
19382 


134  3 

1 944 
J9514 
1 9586 
19682 


20C 


20282 
~ 


32 

33 

34 

35 

36 


37 

38 


SW 


I F( F5-FCAL) 34, 34, 33 
PRINT  72 
GO  TO  1 

IF  ( CH2-SSPAR) 35 , 36,36 
PRINT  73 
GO  TO  1 
PRINT  74 
DO  38  J= 1 , N 
DO  37  1=1,4 

YB( S )=(10.**((C( I ) — A ( J ) )/B ) ) /FAL 
PRINT  63, IDE,J,YB(1),YB(2),YB(3),YB(4) 

PRINT  51 
PRINT  75 
DO  39  J=2,N 
A1=(A(J)-A(1))/B 
Q=A1 -XB ( 1 ) +XB ( J ) 

E=A1+(G*Q/(1 .-G)) 

A3=1 . -G 

A3=((T5*Z)/(B*A3))*SQRT((A3/SWN(1))+(A3/SWN(J))+(Q*Q/SSX)) 
P0TUL=1 0.**(E+A3) 

P0T=1 0.**A1 
P0TLL=1  0.**(E-A3) 

PRINT  63, IDE, J,POTUL,POT,POTLL 
PRINT  51 
GO  TO  1 
END 

OFF  TO  IGNORE  SUBROUTINES 
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10 


Do  13 


11 


\2 


Do  j=l,N 


o\ 

i — ! 

II 

•r-j 

o 

Q 


Z = 1,00 
T05  « 1,96 


gjgg; 


Z = x/SDEV 


STDEB  = VZ2/SSX 
G *=  T05*T05*Z/SREG 


IT,T05,STDE,STDEB,B 


A 


YB(i)  = (IO.**((C(!)-A(j))/B))/FAL  I 


H i 
u 

o (print  y (jDE,j,YB(l),YB(2),YB(3),YB(U)  ) 


■^7 
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POTENCY  PROBIT  ANALYSIS  MAXIMUM  LIKELIHOOD 


16 


i DENT 

SET 

DOSE 

NET  RESPONSE 

: CHECK 

1 01 

1 

4, 000000 

1 6 . 666666 

.000000 

1 01 

1 

8,000000 

40.000000 

.000000 

1 01 

1 

24.000000 

65.000000 

.000000 

1 01 

1 

48.000000 

78.500000 

.000000 

1 01 

2 

4.000000 

2.666666 

.000000 

1 01 

2 

8.000000 

5.714285 

.000000 

1 01 

2 

24.000000 

28 . 444444 

.000000 

1 01 

2 

48.000000 

45.777777 

.000000 

1 01 

3 

24.000000 

30,000000 

.000000 

1 01 

3 

48,000000 

39.000000 

.000000 

SET  NO. 

DOSES 

SSX 

SSY 

SSXY 

SSREG 

1 

b 

63.918185 

154.699630 

98.430363 

151.577140 

2 

b 

32.391504 

1 05.296120 

58.167962 

1 04.45  676  0 

3 

2 

2.649517 

1 .705918 

2.125996 

1.705918 

SET  NO. 

DOSES 

XBAR 

YBAR 

SNW 

INTERCEPT 

1 

b 

1 . 189489 

5.08  0578 

446.358190 

3.172714 

2 

b 

1 .403076 

4.419619 

332.977760 

2.169175 

3 

2 

1.541593 

4.609921 

117.564690 

2.137305 

ANALYSIS  OF 

LINEAR  REGRESS! 

ION 

SOURCE 

VARIATION  DF 

SS 

MS 

TOTAL 

7 

261 .701660 

7.710  F05 

REGRESSION 

1 

254.583740 

254.583740 

257.035  FCAL 

DEV  REGRESS 

ION 

4 

3.961850 

.99  046  2 

9.490  CHI  SQ 

PARALLELISM 

2 

3. 156070 

1 .578  035 

5.990  CHI  SQ 

IT 

G 

T05 

STDE 

STDEB  SLOPE 

4 

.015089 

1.960000 

1.000000 

.100524  1.603936 

IDENT 

SET 

LD30 

LD50 

LD70 

LD90 

1 01 

1 

6.490901 

13.78  0042 

29.254731 

86.753000 

1 01 

2 

27.414292 

58.199950 

123.557240 

366.400940 

1 01 

3 

28.697683 

60.924561 

129.341500 

383.553880 

IDENT 

SET 

UPPER  LIMIT 

POTENCY 

LOWER  LIMIT 

1 01 

2 

.299413 

.236770 

. 187233 

1 01 

3 

.30631  1 

.226182 

.167013 

POTENCY  PROBIT  ANALYSIS  MAXIMUM  LIKELIHOOD 
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DENT 

SET 

DOSE 

NET  RESPONSE 

CHECK 

1 06 

1 

4 . 000000 

34 . 666666 

.000000 

1 06 

1 

8 . 000000 

71 .500000 

.000000 

1 06 

1 

24 . 000000 

86. 153846 

.000000 

1 06 

1 

48  * 000000 

96.500000 

.000000 

1 06 

2 

4.000000 

17.333333 

.000000 

1 06 

2. 

8.000000 

40.571428 

.000000 

106 

2 

24.000000 

77.333333 

.000000 

1 06 

2 

48.000000 

89.333333 

.000000 

1 06 

3 

24.000000 

60.000000 

.000000 

1 G6 

3 

48.000000 

69.000000 

.000000 

SET  NO. 

DOSES 

SSX 

SSY 

SSXY 

SSREG 

1 

4 

45.918937 

165.976230 

84.454855 

155.330740 

2 

4 

45.486714 

182.91151 0 

91  .0178  00 

182. 124380 

3 

2 

2.645250 

1 .613487 

2.065932 

1 .613487 

SET  NO. DOSES 

XBAR 

YBAR 

SNW 

INTERCEPT 

1 

4 

1 .032328 

5.572027 

345.789870 

3.623315 

2 

4 

1.213743 

5.367373 

347.655920 

3.0762  03 

3 

2 

1 .518727 

5.359741 

1 17-51  0200 

2.492857 

ANALYSIS  OF  LINEAR  REGRESSI 

ION 

SOURCE 

VARIATION 

1 DF 

ss 

MS 

TOTAL 

7 

350.501220 

7.710  F05 

REGRESSION 

1 

335.1371 00 

335.137100 

117.256  FCAL 

DEV  REGRESSION 

4 

1 1 .432620 

2.858155 

9.490  CHI  SQ 

PARALLELISM 

2 

3.931500 

1 .965750 

5.990  CHI  SQ 

IT 

G 

TQ5 

STDE 

STDEB  SLOPE 

3 

. 065720 

2.7/6000 

1 .690607 

.174325  1.887686 

IDENT 

SET 

LD30 

LD50 

LD70 

LD90 

1 06 

1 

2.828115 

5.361638 

1 0.164779 

25.599054 

1 06 

2 

5.512247 

1 0.45031  0 

19.812055 

49.894826 

1 06 

3 

11.229392 

21 .28  9 066 

40.3  6 0545 

1 01 .644290 

IDENT 

SET  UPPER  LIMIT 

POTENCY 

LOWER  LIMIT 

1 06 

2 

.794596 

.513060 

.319843 

1 06 

3 

.467314 

.251849 

. 13  08  7 3 

Code 


Definition  of  Code  Used  in  Printed  Output 
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IDEM* 

SET 


DOSE 

NET  RESPONSE 
CHECK 
HO  DOSES 

ssx 


SSY 

SSXY 


SSREG 


XRAK 

YEAR 

snw 

INTERCEPT 

TOTAL 

REGRESSION 


DEV  REGRESSION 


PARALLELISM 


DF 

F05 

FCAL 
CHI  SQ 

CHI  SQ 


Identification  of  the  experiment  or  set  of  data. 

Identification  of  subset  of  data  - subsets  should  be  numbered 
consecutively  for  consistency  of  identifying  numbers . 

The  dose  administered. 

The  percent  response  adjusted  for  response  in  untreated  control. 
(Note  each  subset  may  have  a different  control.) 

The  -proportion  that  responded  in  the  untreated  control , ( Note 

each  subset  may  have  a different  control . ) 

The  number  of  doses  in  each  subset  of  data. 

The  weighted  sum  of  squares  of  log  dose  corrected  for  its  mean  for 
each  subset  of  data. 

The  weighted  sum  of  squares  of  probit  response  corrected  for  its 
mean  for  each  subset  of  data. 

The  sum  of  the  weighted  cross  products  of  log  dose  and  probit 
response  corrected  for  their  means  for  each  subset  of  data. 

The  sum  of  squares  of  probit  response  attributable  to  linear 
regression  of  probits  on  log  dose  for  each  subset  of  data. 

The  weighted  mean  of  log  dose  {£  * x bar) , the  mean  may  be  nega- 
tive if  doses  were  less  than  one  and  ’’log  factor”  not  used. 

The  weighted  mean  of  probit  response. 

The  sum  of  the  weighting  coefficients  for  each  subset  of  data. 

The  intercept  of  each  regression  line,  or  the  value  of  Y (probit 
response)  when  x (log  dose)  is  zero  (0). 

The  sum  of  the  SSY 0 s ( see  above } . 

The  sum  of  SSXY  squared  and  divided  by  the  sum  of  the  SSX  * s , 
which  is  the  sum  of  squares  attributable  to  a single  slope  for 
all  lines  (subsets  cf  data). 

The  sum  of  the  SSI's  minus  the  sum  of  the  SSREG* s (see  above, 
but  do  not  confuse  with  REGRESSION),  which  is  the  weighted  sum 
of  squares  of  the  deviations  from,  parallel  regression  lines 
corrected  for  their  means  (i *e« , the  within  deviations  sum  of 
squares  ) . 

The  weighted  sum  of  squares  attributable  to  deviations  from 

parallelism,  which  is  the  difference  between  the  sum  of  SSREG’s 
and  the  sum  of  5SXY"s  squared  and  divided  by  the  sum  of  SSX’s 
which  is  also  the  difference  in  sum  of  squares  resulting  from 
fitting  one  slope  to  all  subsets  of  data  and  from  fitting  indivi- 
dual slopes  for  each  subset  of  data.  1/ 

Degrees  of  freedom  associated  with  each  source  of  variation. 

The  critical  value  of  F at  5$  level  with  1 and  2 (k^-l)-N  degrees 
of  freedom,  where  is  the  number  of  doses  in  each  of  the  i_ 
subsets  of  data  and  N is  the  number  of  subsets  of  data  or 
number  of  lines. 

The  calculated  F values  from  test  of  significance  of  the  regression 
which  is  the  ratio,  REGRESSION /BIT*?  REGRESSION  mean  squares. 

Chi  square  value  at  5$  probability  level  with  2(Kq-l)-N  degrees  of 
freedom. 

Chi  square  value  at  5$  probability  level  with  N-l  degrees  of  free- 
dom. | 


1/  This  definition  is  not  entirely  correct,  but  is  probably  the  procedure  that 
should  be  used.  The  reader  interested  in  exactly  how  the  "parallelism’’  sum 
of  squares  is  calculated  should  examine  the  accompanying  flow  sheet  or  con- 
sult D,  J.  Finney’s,  Probit  Analysis. 
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Code 


Definition  of  Code  Used  in  Printed  Output 


IT 


G 


TO  5 


STDS 


STDKB 

SLOPE 

LB -values 


The  number  of  iterations  required  to  obtain  a difference  between 
the  last  and  the  next  to  last  estimates  of  the  slope  (B)  of  less 
than  .00005-  The  number  of  iterations  is  counted  starting  from 
the  second  estimate  of  B as  1.  The  first  estimate  of  B is  ob- 
tained from  a weighted  linear  regression  of  empirical  probits 
on  log  dose , in  which  the  net  responses  are  used  in  the  weighting 
coefficients  rather  than  the  proportions  corresponding  to  the 
fitted  provisional  probits. 

The  precision  of  estimating  the  slope  (will  normally  be  the  ratio 
of  F05/PCAL).  If  G is  less  than  .05,  G is  set  equal  to  zero  in 
calculating  confidence  limits . 

The  critical  value  of  students  t at  5$  probability  level  with 
Z (K^-l)-N  degrees  of  freedom  which  is  used  in  setting  confidence 
limits  about  the  potencies.  If  data  are  heterogeneous  T05 
with  (K^-l)-N  degrees  of  freedom  will  appear  here.  If  data  are 
homogeneous  1.96  will  appear  here. 

The  standard  error.  If  the  data  are  homogeneous,  the  value 
1.000000  will  appear  here.  If  the  data  are  heterogeneous, 
the  value  \/BEV  REGRESSION  mean  square  will  appear  here. 

The  standard  error  of  the  slope  = \J  STDE^/  SSX ' s . 

The  common  slope  of  all  regression  lines  for  the  data  composing 
"experiment . " 

The  ID-30,  ID-50,  LD-70,  and  ID-90  for  each  computed  line.  The 
computed  lines  may  be  obtained  by  plotting  these  values  against 
30,  50,  70,  and  90$  values  on  log-probit  paper.  The  observed 
percent  responses  may  also  be  plotted  to  give  a graphical  pre- 
sentation of  the  numerical  results. 

The  upper  95$  confidence  limit  of  potency. 

The  ratio  of  each  succeeding  LD-50  to  the  first  LD-50  listed. 

The  lower  95$  confidence  limit  of  potency. 


UPPER  LIMIT 
POTENCY 
LOWER  LIMIT 
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INSTRUCTIONS  FOR  PREPARING  PARAMETER  CARDS 


One  parameter  card  is  required  for  each  subset  (up  to  five)  of  data,.  An  additional, 
parameter  card  for  the  entire  set  of  data  (experiment)  is  also  required  and  precedes 
the  first  parameter  card, 

PARAMETER  CARD  FOR  THE  EXPERIMENT 


The  parameter  card  for  the  set  (experimental)  precedes  all  other  cards  and  is 


referred  to  as  the  parameter  card  for  the  experiment. 

Columns  1-4  Enter  any  4-digit  number  to  identify  the  experiment. 


5-8  Enter  any  4 -digit  number  to  identify  the  first  subset  of  data 

(this  field  may  be  left  blank  but  is  useful  in  ordering  cards 


). 


9-11  Leave  blank:  - not  used. 

12  Enter  number,  N,  of  subsets  of  data  in  the  set  (experiment).  This 
number  instructs  the  computer  to  read  N subsets  of  data. 

13-18  Enter  Chi  square  value  at  Jp  probability  level  with  2 (K. -l)-N 


degrees  of  freedom,  where  is  the  number  of  doses  in  each  of 
the  i subsets  of  data,  and  N the  number  of  subsets  in  the  set 
(experiment).  The  Chi  square  value  may  be  entered  any  place 
within  this  six  digit  field  but  decimal  must  be  punched  (see 


NOTE  at  end  of  instructions ) . 

19-24  Enter  Chi  square  value  at  5 °h  probability  level  with  N-l  degrees 
of  freedom.  The  Chi  square  value  may  be  entered  any  place 
within  this  six  digit  field  but  decimal  must  be  punched  (see 


NOTE  at  end  of  instructions ) . 

25-30  Enter  critical  value  of  Student’s  t with  2(K.;-l)-N  degrees  of 
freedom  at  yjo  probability  level.  Decimal  must  be  punched. 

31-3o  Enter  critical  value  of  F at  5$  probability  level  with  one  and 
2(Kq-l)-i<  degrees  of  freedom.  Decimal  point  must  be  punched. 

37-42  Enter  LOG  FACTOR  - 10^  where  i_  is  an  integer  between  0 and  6 

such  that  when  the  lowest  dose  is  multiplied  by  10  ^ it  will  be 
greater  than  unity.  If  all  doses  are  already  greater  than 
unity  the  LOG  FACTOR  will  be  10c= 000001,  If  the  lowest  dose 
is  .00005  the  LOG  FACTOR  will  be  10 5 = 100000,  This  LOG 
FACTOR  eliminates  the  possibility  of  a negative  mean  for  log 
dose,  which  may  confuse  the  biologist.  Negative  values  will 
not  effect  the  results  noi  interfere  with  the  calculations. 


PARAMETER  CARD  FOR  EACH  SUBSET 


One  parameter  card  precedes  each  subset  of  data. 


Columns  1-4 


5-8 


9-10 

11-12 


Enter  4-digit  number  identifying  the  experiment  (see  above  Column 
1-4). 

Enter  4-digit  number  identifying  the  subset  of  data.  For  con- 
sistency of  output  these  subsets  should  be  numbered  consecutively 
starting  with  1 ( one ) , 

Leave  blank  - not  used. 

Enter  number,  M,  of  doses  or  data  cards  (up  to  twenty)  that  com- 
prise this  subset  of  data.  This  number,  M,  instructs  the  com- 
puter to  read  the  next  M cards  as  data  cards. 
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Columns  13-18  Enter  number  of  animals  observed  in  the  untreated  control  (check) 
for  this  subset  of  data.  Note  that  each  subset  may  have  a 
different  value  for  the  untreated  control  as  well  as  a different 
number  of  doses. 

19-24  Enter  the  number  of  animals  that  responded  in  the  untreated  control 
for  this  subset  of  data.  Columns  13-24  may  be  left  blank  if  no 
check  was  used  or  if  no  animals  responded  in  the  untreated  control. 


POTENCY  PROBIT  ANALYSIS 


INSTRUCTIONS  FOR  PREPARING  DATA  CARDS 

Enter  data  for  each  dose  (dose  used,  number  of  animals  treated,  number  of  animals 
that  responded)  in  the  following  card  columns  using  one  card  for  each  dose. 


Columns 


1-4  Enter  any  4-digit  number  to  identify  the  experiment. 

5-8  Enter  any  4-digit  number  to  identify  the  subset  to  which  these 

data  belong.  For  consistency  in  identification  of  "SET,"  the 
subsets  should  be  numbered  consecutively  beginning  with  1 (one). 

9-14  Enter  dose  or  concentration  used.  The  dose  may  be  entered  any 

place  within  this  field  providing  that  decimal  point  is  punched 
in  its  proper  place.  If  no  decimal  point  is  punched,  the  com- 
puter will  automatically  place  a decimal  point  between  card 
columns  l4  and  15 . 

15-20  Enter  number  of  animals  which  received  the  dose  listed.  If  no 
decimal  point  is  punched,  the  program  will  automatically  place 
a decimal  point  between  card  columns  20  and  21. 

21-26  Enter  the  number  of  animals  that  responded  to  the  dose  listed 
on  this  card.  If  no  decimal  point  is  punched,  the  computer 
will  automatically  place  a decimal  point  between  card  columns 
26  and  27* 


Assemble  all  data  cards  for  each  subset  in  either  ascending  or  decending  order  of 
dose  and  precede  each  subset  with  its  proper  parameter  card  (see  instructions  for 
preparing  parameter  cards ) . 

NOTEs  If  decimal  points  appear  any  place  except  between  the  designated  fields, 
the  decimal  points  should  be  punched.  A punched  decimal  point,  with 
FORMAT  FORTRAN,  will  over-ride  a FORMAT  decimal  point.  Decimal  points 
should  not  be  punched  for  the  numbers  identifying  the  experiment  and  sub- 
sets of  data,  which  are  read  as  fixed  rather  than  floating  numbers  (as 
l4  rather  than  F6.0). 
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POTENCY  PROSIT  ANALYSIS 
Operating  Instructions  for  IBM  1620 


1 . Clear  Memory . 

2.  Load  POTENCY  PROBIT  ANALYSIS  program  and  FORMAT  FORTRAN  Subroutines. 

3.  Load  data  cards  in  following  orders 

a.  Parameter  card  for  '’experiment"  or  for  all  subsets  of  data. 

b.  Parameter  card  for  first  subset  of  data. 

c.  Data  cards  in  ascending  or  decending  order  of  dose  for  first 

subset . 

d.  Parameter  card  for  second  subset  of  data. 

e.  Data  cards  in  ascending  or  decending  order  of  dose  for  second 

subset . 

Additional  subsets,  up  to  five,  may  follow  the  first  parameter  card 
as  outlined  above. 

4.  Additional  sets  (experiments)  may  follow  the  first  as  indicated  in  (3)  above. 

5.  Follow  the  last  subset  with  2 blank  cards. 

NOTES: 

(1)  Sense  switches  are  not  interrogated. 

(2)  When  net  response  is  negative,  the  program  types  out  the  doses  for 

all  subsequent  subsets  within  this  set,  then  goes  to  the  next  set 
of  data  ignoring  the  computations  of  the  rest  of  the  analysis. 

The  decision  to  set  the  check  for  this  subset  in  question  to  zero 
(0)  and  re-run  the  data  is  left  to  the  experimenter  or  biometrician. 
It  may  be  prudent  for  the  biometrician  to  examine  the  reason(s) 
for  different  checks  in  each  subset  before  submitting  the  data, 
and  also  to  examine  the  numbers  used  in  the  checks  for  precision 
of  estimating  the  proportion  that  responded  in  the  untreated  control. 
Abbott's  formula  assumes  that  the  response  in  the  untreated  control 
is  known  without  error. 

(3)  If  a print-out  "LINES  NOT  PARALLEL"  or  "NON  SIG  REGRESSION"  occurs 

the  remainder  of  the  calculations  are  bypassed  and  the  program 
precedes  to  the  next  set  of  data  (experiment),  ignoring  the  remaining 
calculations . 

{h)  Error  F8  in  the  print-out  may  occur  when  the  values  are  too  large 

for  the  space  allotted.  When  this  occurs,  the  parameter  and  data 
cards  should  be  examined  for  order  and  correctness  of  entry  of  the 
values . 
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POTENCY  PROBIT  ANALYSIS 
Data  Used  in  Example  of  Printed  Output 


Number 

Number 

Code 

Set 

Dose* 

Observed 

Responded 

101 

1 

4 

150 

25 

8 

200 

80 

24 

260 

169 

48 

200 

157 

101 

2 

4 

75 

2 

8 

175 

10 

24 

22p 

64 

48 

225 

103 

101 

3 

24 

100 

30 

48 

100 

39 

Check 

— — 

000 

106 

1 

4 

150 

52 

8 

200 

143 

24 

260 

224 

48 

200 

193 

106 

2 

4 

75 

13 

8 

175 

71 

24 

225 

174 

48 

225 

201 

106 

3 

24 

100 

60 

48 

100 

69 

Check 

— — 

000 

*Data  supplied  by  courtesy  of  Dr.  Harrie  M.  Taft  from  a time -mortality  study. 
Boll  weevils  were  exposed  for  4,  8,  24 , and  48  hours  to  treated  foliage.  Set 
one  was  from  "immediately  after  application, " sets  two  and  three  were  from  24 
and  48  hour  weathering  under  artificial  conditions.  Purpose  of  study  was  to 
determine  which  of  the  25  treatments  (all  recommended  for  boll  weevil  control) 
was  the  most  toxic  and  which  lost  its  effectiveness  most  rapidly.  Loss  in 
toxicity  or  loss  in  effectiveness  is  therefore  1. -potency,  and  its  confidence 
limits  may  be  similarly  expressed.  Treatment  101  lost  1.-. 23677  = 76$  of  its 
effectiveness  in  24  hours  and  77$  in  48  hours.  Treatment  106  lost  1.-. 513^6  = 
49$  of  its  effectiveness  in  24  hours  and  75$  in  48  hours.  Treatment  106  was 
originally  the  most  toxic  (LD-50  = 5*36  hours)  while  treatment  101  the  least 
toxic  (LD-50  = 13‘78  hours). 


